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Abstract 


Large  amplitude  wall-pressure  events,  observed  beneath  a  turbulent  bound¬ 
ary  layer,  appear  to  be  the  signatures  of  intermittent  organized  motions 
within  the  turbulent  flow.  The  temporal  localization  of  these  events  could  be 
applied  to  the  active  control  of  turbulent  wall  flows.  This  report  presents  pre¬ 
liminary  results  on  utilizing  time-frequency  localization  techniques  (wavelet 
transforms)  for  the  detection  of  these  events.  The  advantage  of  these  meth¬ 
ods  is  that  they  do  not  require  a  priori  assumptions  regarding  the  features 
of  the  signal.  A  tutorial  overview  of  these  techniques  is  first  presented.  This 
is  followed  by  a  discussion  of  some  exploratory  results  obtained  from  the  ap¬ 
plication  of  wavelet  filtering  to  wall  pressure  and  turbulent  temporal  records 
acquired  from  wind  tunnel  experiments. 


Chapter  1 

INTRODUCTION 


1.1  Background 

It  has  been  widely  postulated  that  large  amplitude  wall-pressure  events  can 
be  used  as  a  non-intrusive  observation  of  the  passage  of  organized  flow  struc¬ 
tures.  These  organized  motions  are  proposed  to  be  the  primary  mechanism 
for  the  production  of  turbulence  (cf.  (!]).  It  is  speculated  that  the  temporal 
localization  of  these  surface  events  could  be  used  for  the  active  control  of 
turbulent  wall  flows. 

The  pressure  fluctuations  beneath  a  turbulent  boundary  layer  are  the  inte¬ 
gral  effects  of  active  (turbulent  producing)  and  passive  (high  kinetic  energy) 
flow  structures  throughout  the  boundary  la3rer.  It  is  the  active  structures 
that  are  of  primary  interest  in  a  control  system.  Therefore,  signal  detection 
techniques  on  the  temporal  records  of  the  pressure  are  required  to  discrimi¬ 
nate  between  the  two  contributions.  In  the  context  of  this  paper,  the  pressure 
fluctuations  of  the  passive  motions  will  be  considered  as  noise  contamination 
on  the  pressure  signatures  of  the  active  motions. 

Previous  work  by  Farabee  and  Casarella  (cf.  [2])  had  shown  that  turbulent 
sources  of  the  wall  pressure  spectrum  can  be  attributed  to  flow  activities  at 
distinct  locations  aEcross  the  boundary  layer.  They  concluded  that  the  inner 
layer  is  characterized  by  high  fi^uency  turbulent  motions  while  the  outer 
layer  is  characterized  by  low  frequency  structures.  It  is  widely  believed  that 
the  near-wall  (active)  structures  can  spatially  extend  to  the  overlap  region 
and  are  thus  not  limited  to  a  distinct  spatial  location.  Clearly,  conventional 
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frequency  filtering  techniques  cannot  be  applied. 

A  more  tractable  approach  is  to  examine  the  intermittent  large-amplitude 
pressure  events  at  the  wall.  These  appear  to  be  the  footprint  of  the  bursting 
process  associated  with  the  near-wall  turbulent  production  mechanism  (cf. 
[3]).  Extensive  work  has  been  done  using  a  range  of  schemes  for  the  simulta¬ 
neous  detection  of  intermittent  large-amplitude  turbulent  and  wall-pressure 
events  (cf.  [4,  5]).  Conditional  sampling  techniques  are  often  used  and  these 
include  detection  algorithms  based  on: 

•  VITA  on  u  with  slope  criteria 

•  VITA  on  u+LEVEL  (cf.  [6]) 

•  Q2  and  Q4  quadrant 

•  Peak  pressure  events 

By  averaging  the  collection  of  detected  events,  a  correlation  can  be  ob¬ 
served  between  turbulent  events  and  peak  pressure  events.  The  number  of 
events  detected,  and  thus  averaging,  are  affected  by  filtering  techniques  and 
various  refinements  to  the  detection  algorithm.  Furthermore,  the  shape  of 
the  detected  events,  including  non-symmetric  patterns,  is  strongly  dependent 
on  the  alignment  procedure  used  for  averaging.  It  appears  that  an  unbiased 
detection  algorithm  is  required.  This  was  the  motivation  for  this  investiga¬ 
tion. 

1.2  Review  of  Temporal  Records 

It  is  necessary  to  first  examine  the  characteristic  features  of  the  temporal 
records  of  both,  the  fluctuating  velocities  within  the  flow,  and  simultaneous 
wall  pressure.  The  data  acquisition  and  conditioning,  spectral  content,  and 
cross-correlation  will  be  presented. 

The  experimental  setup  for  the  measurement  of  the  time  records  used  in 
the  present  study  included  pressure  and  velocity  sensors.  Two  components 
of  the  velocity  (streamwise  and  wall  normal)  were  measured  using  a  cross¬ 
wire  anemometer.  Simultaneous  with  the  velocity  measurements,  the  wall 
pressure  fluctuations  were  recorded  using  a  1/8”  B&K  microphone  with  a 
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Figure  1.1:  The  PSD  of  the  u,  v  and  p  time  records. 

1  /32”  diameter  pinhole  cap.  The  cross-wire  had  a  non-dimensional  sensing 
length  of  22  viscous  units  whereas  the  microphone  had  a  spatial  resolution 
of  34  viscous  units. 

The  presented  data  were  obtained  with  the  hot-wire  located  at  a  distance 
of  y/S  =  0.08  (y"*"  =  84),  flow  speed  of  51  ft/s  (15.5  m/s)  and  =  2945 
(<$^  =  1221).  A  detailed  description  of  the  wind  tunnel  facility  and  the 
measuring  procedures  can  be  found  in  [7]. 

The  time  records  of  the  wall  pressure  (p),  the  streamwise  velocity  (u)  and 
the  wall-normal  velocity  (v)  were  obtained  by  antialiasing  the  sensor  signals 
with  8-pole  low-pass  elliptic  Alters  with  cutoff  frequency  at  12.5  KHz.  The 
data  were  digitized  at  a  sampling  rate  of  25  KHz  for  10  seconds.  These 
data  records  were  then  filtered  and  normalized.  The  filtering  was  performed 
digitally  using  a  5th  order  Butterworth  high-pass  filter  with  the  cutoff  fre¬ 
quency  at  100  Hz.  Wilczynski  (cf.  [7])  showed  that  high-pass  filtering  at  this 
frequency  did  not  affect  the  integrity  of  the  large  amplitude  events.  The  nor¬ 
malization  was  performed  with  respect  to  the  power  present  in  the  resulting 
sequences,  so  that  the  three  normalized  time  records  have  unit  power. 

The  estimates  of  the  power  spectral  density  (PSD)  for  p,  u  and  uv  are 
given  in  figures  1.1  and  1.2.  These  estimates  were  computed  using  Welch’s 
averaged  periodogram  method.  The  frequency  bump  in  the  p  spectrum  is 
due  to  the  shedding  from  the  hot-wire  traverse  system. 

The  higher  frequency  components  of  the  uv  record  are  greater  than  their 
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Figure  1.2:  The  PSD  of  the  u,  v  and  uv  time  records. 
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Figure  1.3;  The  biased  cross-correlations  between  u  and  v,  u  and  p,  and  p 
and  V.  These  computations  were  performed  with  32768  samples  of  each  time 
record. 
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equivalent  components  in  the  u  and  v  records.  This  is  a  consequence  of  the 
increase  of  axis  crossings  in  uv.  Note  that  the  roll-off  followed  by  the  PSD’s 
of  the  time  records  are  similar. 

The  physics  of  the  turbulent  flow  can  be  initially  seen  in  the  cross- 
correlation,  i2iy  =  E{x{ti)y{t2)},  between  the  u,  v  and  p  time  records,  as 
shown  in  figure  1.3.  The  biased  estimates  of  the  cross-correlations  plotted  in 
figure  1.3  show  that  u  is  correlated  to  —v  and  p  only  at  ti  —  t2  ~  0.  It  also 
shows  that  p  and  v  are  correlated.  The  cross-correlation  results  are  given  as 
a  function  of  =  tul/v  where  t  is  the  sampling  time.  A  more  extensive 
discussion  of  the  correlation  based  on  cross-spectra  data  can  be  found  in  [5] 
and  [7]. 


1.3  Objective  of  the  Investigation 

A  formulation  of  the  problem  will  be  presented  prior  to  the  application  of 
time-frequency  localization  techniques. 

The  u,  V  and  p  time  records  can  be  assumed  to  be  random  processes  of 
the  form 


u(<) 

=  «,(<) -f- u„(t) 

(1.1) 

u(t) 

=  Vs{t)  +  Vnit) 

(1.2) 

pit) 

=  Psit)+Pnit) 

(1.3) 

where  u,(t),  u,(t)  and  pj(t)  are  due  to  active  structures  in  the  turbulent 
flow,  and  u„(t),  and  p„(<)  are  assumed  to  be  noise  (in  general,  non- 
gaussian).  It  is  postulated  that  u,  and  w,  represent  the  distinct  features  of 
organized  structures  and  that  ps{t)  represents  the  unique  signatures  of  these 
intermittent  flow  events. 

Different  schemes  could  be  exercised  to  isolate  the  signal  from  the  noise 
in  the  time  records.  This  study  uses  time-frequency  localization  procedures 
in  an  attempt  to  partition  the  data  as  described  by  equations  1.1  thru  1.3. 

Chapter  2  will  examine  the  application  of  three  techniques  for  the  char- 
aicterization  of  time  records;  Short-time  Fourier  Transforms,  Continuous 
Wavelet  Transforms,  and  Discrete  Orthogonal  Wavelet  Transforms.  A  filter¬ 
ing  technique  for  the  extraction  of  the  signal,  based  on  an  orthogonal  wavelet 
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expansion  of  the  temporal  records  is  presented  in  Chapter  3.  A  comparative 
evaluation  of  event  detection  between  the  filtered  and  unfiltered  temporal 
records  is  briefly  presented  in  Chapter  4.  Finally,  a  discussion  and  summary 
of  the  findings  are  included. 
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Chapter  2 

TIME-FREQUENCY 

LOCALIZATION 


The  segments  of  the  time  records  shown  in  figure  2.1  illustrate  the  diffi¬ 
culty  of  objectively  identifying  organized  structures  that  are  common  to  all 
the  time  records.  It  is  possible  that  in  transforming  the  time  records  to 
the  frequency  domain,  or  to  the  time-frequency  domain,  the  structures  will 
become  more  clearly  identifiable.  There  are  several  techniques  for  implement¬ 
ing  these  Jilternative  representations.  The  Short  Time  Fourier  Transform  and 
the  Wavelet  Transform  will  be  examined  in  this  paper.  These  techniques  will 
be  applied  to  the  u{t)  record.  Similar  results  have  been  obtained  for  v{t), 
uv{t)  =  u{t)v{t)  and  p{t)  records. 


2.1  Short  Time  Fourier  Transform 

The  Short  Time  Fourier  Transform  (STFT)  (cf.  [8])  defined  as 

STFT(t,  f)  =  J  /(r)j(T  -  i)t-‘^>’dT  (2.1) 

(with  g{t)  assumed  as  a  Hanning  window)  is  applied  to  the  interval  of  u  shown 
in  figure  2.1.  The  magnitude  and  phase  spectrograms  of  the  corresponding 
STFT  are  shown  in  figure  2.2. 

Clearly  identifiable  structures  are  absent  in  the  spectrograms  of  figure  2.2. 
This  is  also  the  case  for  other  sections  randomly  selected  from  the  time 
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records.  The  magnitude  spectrogram  shows  that  the  main  frequency  compo¬ 
nents  in  the  time  record  are  located  below  4  KHz,  which  is  already  known 
from  the  PSD  plots  in  figures  1.1  and  1.2. 

Windowed  spectra  of  conditionally  sampled  events  have  also  been  tried 
without  any  new  insight  on  the  detection  problem.  The  results  only  con¬ 
firm  the  broad  spectral  range  of  short-time  events.  Prom  these  preliminary 
studies,  STFT  techniques  do  not  appear  to  be  a  feasible  means  for  signal 
identification. 


2.2  Continuous  Wavelet  Transform 

A  comprehensive  review  on  the  theoretical  foundation  of  wavelets  is  given 
by  Daubechies  (cf.  [8]),  and  formed  the  basis  of  much  of  the  work  to  be 
presented.  The  applications  of  wavelet  transforms  to  fluid  mechanics  and 
turbulence  had  been  exantined  by  Meneveau  [9]  ,Farge  [10]  and  Zubair  [11]. 

The  Continuous  Wavelet  Transform  (WT)  is  defined  as 

WT(a,6)  =  J  S{T)r'\r)dr  (2.2) 

where  the  are  known  as  wavelets  and  xl){t)  is  known  as 

the  mother  wavelet  since  all  waveletes  are  mutually  similar  with  a  dilation  or 
scaling  parameter  o.  Individual  wavelets  are  characterized  by  three  features; 
shape,  duration  and  location  (cf.  [12]).  A  wide  range  of  both  continuous  and 
discrete  wavelets  exists  and  the  selection  dictates  the  shape.  The  parameter 
6  defines  the  location  along  the  temporal  axis,  while  o  is  associated  with 
the  duration.  It  should  be  noted  that  decreasing  values  of  a  correspond  to 
wavelets  with  higher  amplitude  and  shorter  duration. 

An  equivalent  expression  for  the  Wavelet  Transform  can  be  formulated  in 
terms  of  the  Fourier  TVansforms  for  the  signal  (/(u;))  and  wavelets  (^(ow)) 
as  follows 

WT{a,b)  =  y/a  f  f{uj)'0{au))e^^duj.  (2.3) 

J-OO 

This  expression  illustrates  that  the  wavelet  decomposition  acts  as  a  filtering 
of  f{uj)  by  v^^(aa))  with  a  phase  lag  ftw.  For  most  practical  applications, 
wavelets  are  well-localized  in  both  time  and  frequency.  In  addition,  wavelet 
filters  have  constant  relative  bandwith  {Au/uj  =  Constant)  in  contrast  to 
the  constant  bandwidth  of  the  STFT. 
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A  discussion  comparing  STFT  and  Wavelet  Transforms  is  presented  by 
Daubechies  (cf.  [8]).  She  explains  that  the  essential  difference  between  wavelets 
and  windowed  Fourier  transforms  is  in  the  shape  of  the  windowed  function. 
The  STFT  consists  of  a  g{t)  function  with  constant  bandwidth,  translated  to 
the  proper  location,  and  filled  in  with  high  frequency  oscillations.  Wavelets 
are  a  family  of  prescribed  shapes  that  are  also  translated,  but  com¬ 
pressed  or  stretched.  The  time  widths  are  adapted  to  their  frequency  bands; 
wavelets  with  high-frequency  components  are  very  narrow  while,  wavelets 
with  low-frequency  components  are  much  broader.  As  a  result,  the  wavelet 
transform  is  better  able  than  the  windowed  Fourier  transform  to  zoom  in  on 
very  short-lived  high  frequency  phenomena,  such  as  transients  in  signals. 

The  scalograms  shown  in  figure  2.3  were  computed  from  the  u(i)  data 
record  using  (2.2)  and  the  Morlet  wavelet  with 

a  =  6/5,  C  =  1  and  A  :§>  2/ir.  These  scalograms  do  not  show  distinct  concen¬ 
trations  of  energy  that  could  be  unequivocally  identified.  Real  wavelets,  like 
the  First  Derivative  of  the  Gaussian,  ipit)  =  — or  the  Second  Deriva¬ 
tive  of  the  Gaussian  (known  also  as  the  Mexican  Hat),  tj;  =  {1  — 
give  equivalent  results  as  illustrated  in  figure  2.4.  Liandrat  et  al  [13]  obtained 
similar  results  to  those  shown  in  figure  2.4  for  the  Mexican  Hat.  They  sug¬ 
gested  that  shaded  cone-like  structures  in  the  scalogram  correspond  to  ejec¬ 
tions/sweep  events.  A  comparison  between  their  wavelet  results  and  VITA 
detection  showed  reasonable  agreement. 


2.3  Discrete  Orthogonal  Wavelet  Transform 

2.3.1  Basic  Concepts 

Similar  to  the  Fast  Fourier  Transform  (FFT),  Discrete  Orthogonal  Wavelet 
Transform  can  be  formulated.  For  this  case,  a  particular  set  of  wavelets  is 
specified  by  a  discrete  set  of  wavelet  (filter)  coefficients.  Several  families  of 
orthogonaJ  wavelets  with  compact  support  have  been  defined  in  the  literature 
Fast  algorithms  have  been  developed  for  the  computation  of  these  wavelets 
(cf.  [14]).  Zubair  [11]  utilized  a  new  class  of  algorithms,  called  wavelet- 
packets.  This  procedure  segments  the  turbulence  data  adaptively  and  selects 
the  best  basis  from  a  large  collection  of  wavelets.  Intuitively,  it  represents 
the  wavelets  expansion  with  the  least  number  of  coefficients. 
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Figure  2.3:  Scalograms  of  the  magnitude  and  phase  of  the  Morlet  WT  of  u{t) 
for  2093  <  <  2360.  The  Morlet  wavelet  has  a  =  6/5  and  C  =  1.  The 

maximum  and  minimum  intensity  values  for  the  magnitude  plot  are  2.4e+05 
and  32.  The  maximum  and  minimum  intensity  values  for  the  phase  plot 
are  tt  and  — tt. 

This  study  is  restricted  to  the  application  of  three  wavelet  families,  two 
of  them  constructed  by  I.  Daubechies  {Extremal  Phase,  and  Least  Asymmet¬ 
ric),  and  one  formulated  by  R.  Coifman  and  documented  by  I.  Daubechies 
{Coiflets)  (cf.  [8]).  These  are  denoted  as: 

•  D8  —  Daubechies  Extremal  Phase  8-coefficients, 

•  L8  —  Daubechies  Least  Asymmetric  8-coefficients, 

•  Cl8  —  Coiflets  18-coefficients. 

The  Fast  Wavelet  Transform  (FWT)  of  a  one  dimensional  time  record 
containing  2^  samples  is  a  vector  of  length  2^.  In  order  to  compute  a  FWT, 
the  number  of  samples  in  the  time  record  must  be  a  power  of  two  (cf.  [14]). 
This  is  not  a  severe  limitation  because  the  time  records  can  be  conveniently 
zero  padded.  The  u,  v  and  p  records  used  in  this  analysis  contain  249000 
samples  (digitized  at  25000  samples  per  second).  The  corresponding  zero 
padded  sequences  contain  2^®  =  262144  samples. 
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Figure  2.4:  Scalograms  of  the  First  and  Second  Derivatives  of  the  Gaussian 
WT  of  u(t)  for  2093  <  <  2360.  The  maximum  and  minimum  intensity 

values  for  the  First  Derivative  WT  plot  are  426  and  -588.  The  maximum 
and  minimum  intensity  values  for  the  Second  Derivative  WT  plot  are  335 
and  -290. 

Discrete  Orthonormal  wavelets  with  indices  m,  n  are  defined  as 

r'-’it)  =  =  y^V'(2-"(  -  n2-”‘)  (2.4) 

where  Om  =  2"*  and  6„  =  n. 

A  function  f{t)  is  represented  by  the  wavelet  expansion  with  coefficients 
Ct  where  0  <  £  <  2^  —  1.  The  index  i  is  related  to  m,  n  by  means  of  a 
dyadic  pattern  in  the  m,  n  plane  (figure  2.5).  There  are  TV  rows  of  coefficients 
denoted  by  1  <  m  <  N.  For  each  row  m  there  are  2^~”*  coefficients,  therefore 
the  range  of  n  in  that  row  is  1  <  n  <  2^“’”.  The  index  £  =  1  corresponds 
to  m  =  TV  and  n  =  1,  the  single  element  row.  When  displaying  the  wavelet 
coefficients  with  increasing  index  t,  they  are  grouped  in  row  intervals  2^^“^^  < 
£  <  2*  —  1  with  k  =  N  ~  m  +  \  and  1  <  A:  <  TV.  The  number  of  terms  in 
each  row  k  doubles  sequentially. 
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Figure  2.5:  Dyadic  Arrangement  of  Wavelet  Coefficients 

2.3.2  Example  using  the  D8  wavelet  transform 

The  features  of  discrete  wavelets  and  the  dyadic  grid  expansion  will  be  il¬ 
lustrated  primarily  for  the  D8  wavelets  with  a  data  record  of  length  2^  for 
iV  =  18. 

Figure  2.6  shows  some  of  the  orthogonal  wavelet  functions  for  different 
values  of  m  and  n.  It  should  be  observed  that  the  scaling  factor  Om  =  2”* 
controls  not  only  the  time  duration  of  the  wavelet,  but  also  its  energy.  There¬ 
fore,  the  peak  amplitude  of  short  wavelets  (for  instance,  m  =  4  )  is  higher 
than  the  peak  amplitude  of  long  wavelets  (for  instance  m  =  6).  The  normal¬ 
ized  magnitudes  of  the  FFT’s  corresponding  to  the  wavelets  in  figure  2.6  are 
shown  in  figure  2.7  and  illustrates  the  constant  relative  bandwidth  (Aa’/u;) 
of  the  wavelet  filters.  It  should  also  be  noted  that  the  sequential  terms  in  the 
wavelet  expansion  will  have  overlapping  frequency  bands,  particularly  in  the 
high  frequency  groups.  This  is  in  contrast  to  traditional  band-pass  filtering 
techniques  [14]. 

For  this  example  (with  N  =  18),  the  dyadic  arrangement  of  the  coeffi¬ 
cients  for  1  <  f  <  262143  can  be  sununarized  in  tabular  form.  Table  2.1 
lists  the  number  of  terms  in  each  of  the  m  rows  of  the  dyadic  pattern;  the 
center  frequency  fc  (computed  for  D8  wavelets)  for  those  terms  in  the  respec¬ 
tive  rows;  and  the  time  interval  T'*'  between  the  n  wavelets  stretched  along 
the  time  record  (as  shown  in  figure  2.6).  It  should  be  noted  that  the  scaling 


Figure  2.6:  Daubechies  Extremal  Phase  orthogonal  base  functions  computed 
with  FIR  filters  having  8  coefficients. 


Figure  2.7:  Normalized  FFT  mi^itude  of  the  Daubechies  Extremal  Phase 
orthogonal  base  functions  of  figure  2.6. 
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Table  2.1;  Grouping  of  wavelet  coefficients  for  a  dyadic  grid  with  N  =  18. 
The  energy  of  u{t),  v{t)  and  p{t)  records  is  concentrated  in  bands  marked 
with  •  symbol 
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factor  a„  =  2"*  is  equal  to  the  time  interval  T* . 

Let  be  the  Fourier  Transform  of  the  mother  wavelet.  Then  it  can 
be  shown  that  the  Fourier  IVansform  of  is  proportional  to 

It  is  clear  that  the  center  frequency  fc  of  a  wavelet  in  row  m  is  related  to 
the  center  frequency  /o  of  the  mother  wavelet  through  fc  =  /o/2'''  =  JolT*  ■ 
Since  m  =  N  -  k  +  1,  then  fc  =  (2*/2^'*’‘)/o  or  fjf,  =  >12*  where  A 
is  a  proportionality  constant  and  /,  is  the  sampling  frequency.  Thus  fc 
doubles  sequentially  for  each  row  k  as  shown  in  Table  2.1.  This  sequence  of 
frequencies  can,  to  a  first  approximation,  be  equated  to  the  Fourier  frequency. 
The  assumption  being  that  the  overlapping  bands  are  not  significant. 

As  an  illustration  of  the  magnitude  of  the  coefficients  for  a  typical  signal, 
the  u(t)  temporal  record  was  analyzed  by  means  of  the  wavelet  expansion. 
Two  wavelet  families  were  examined;  Dfi-Daubechies  Extremal  Phase  and 
L8-Daubecies  Least  Asymmetric.  For  purposes  of  graphical  display  =  15 
was  chosen.  Figure  2.8  shows  the  time  record  and  the  magnitudes  of  the 
coefficients  for  1  <  ^  <  32767  for  D8  and  L8  wavelet  expansions.  The  last 
row  with  the  highest  concentration  of  coefficients  (m  =  1  or  /:  =  15)  has  the 
coefficient  indices  2**  <  ^  <  2^®  - 1.  The  magnitudes  of  these  coefficients  are 
relatively  small. 

It  will  be  shown  in  subsequent  results  that  only  approximately  10%  of 
the  coefficients  are  needed  to  represent  the  signal  and  for  this  application 
are  dominated  by  terms  in  the  expansion  for  2  <  m  <  8.  This  covers  the 
center  frequency  range  34  <  /c  <  2181,  which  is  not  surprising  in  view  of  the 
spectral  content  of  the  time  records  shown  in  figure  1.1. 

The  energy  of  the  time  record  can  be  computed  from  the  wavelet  coeffi¬ 
cients  by  application  of  Parseval’s  theorem.  When  the  basis  wavelet  functions 
are  orthonormal,  the  expected  value  of  the  signal  energy  is  given  by  the  mean 
square  value  of  the  wavelet  coefficients.  The  Power  contauned  in  each  row 
of  the  expansion  can  be  correlated  with  the  center  frequency  of  the  wavelets 
contained  in  the  respective  rows.  Figure  2.9  displays  the  Power  versus  fc 
for  the  D8  wavelets.  As  shown  in  the  figure,  the  shape  of  this  curve  does 
not  approximate  the  spectral  curves  of  u(t).  These  results  offer  a  poor  ap¬ 
proximation  to  the  spectral  content  of  the  time  records  because  the  center 
frequency  only  represents  the  location  of  the  frequency  bands  which,  as  pre¬ 
viously  discussed,  have  extensive  overlapping. 

A  more  meaningful  wavelet  spectrum  that  accounts  for  the  overlap  has 
been  formulated  by  Zubair  [11].  The  energy  spectrum  of  the  wavelet-coefficients 
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Figure  2.9;  Wavelet  Power  versus  Center  Frequency  for  the  D8  wavelets  and 
the  u{t)  record. 
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is  segmented  into  bins.  The  wavelet-spectrum  is  found  by  adding  the  energy 
contributions  from  all  the  wavelet  coefficients  contained  in  each  of  the  bins. 
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Chapter  3 

APPLICATION  OF 
WAVELET  FILTERING 


Traditional  bandpass  filtering  techniques  based  on  Fourier  decomposition 
have  often  been  used  to  investigate  intermittent  turbulent  structures.  These 
techiques  have  had  limited  success  and  often  present  a  biased  viewpoint  on 
the  scales  of  turbulence.  Wavelet  filtering  offers  a  promising  new  approach 
to  this  problem. 

Zubair  [11]  examined  the  application  of  wavelet  transforms,  including 
spectral  analysis  and  filtering,  to  several  aspects  of  turbulent  flows.  He  de¬ 
veloped  a  wavelet  power-spectra  and  filtering  technique  analogous  to  Fourier 
techniques  to  interpret  the  scale-based  behavior  of  turbulence.  His  conclu¬ 
sion  is  that  wavelet-packet  filtering  is  an  improvement  over  Fourier  filtering 
by  providing  sharper  separation  between  bands  of  scales,  and  thus  offering 
more  sensitivity  than  Fourier  filtering.  However,  his  results  on  three  facets 
of  the  structure  of  turbulence  —  local  itotropy,  intermittency  and  scaling  of 
structure  functions  —  showed  only  modest  improvement  over  Fourier  filter¬ 
ing. 


3.1  Conventional  Wavelet  Filtering 

A  wavelet  filter,  based  on  an  orthogonal  wavelet  expansion,  can  be  imple¬ 
mented  by  first  sorting  the  wavelet  coefiScients  in  decreasing  magnitude  order; 
and  then  zeroing  out  all  the  coefficients  with  magnitude  lower  than  a  speci- 
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fied  threshold  (cf.  [14]).  The  limit  can  be  set  so  that  the  energy  contained  in 
the  filtered  sequence  is  a  fraction  of  the  total  energy.  Let  kj  he  the  fraction 
of  the  total  signal  energy  that  is  present  in  the  filtered  sequence,  then 

Eivlit)}  =  kfE{y\t)}  (3.1) 

where  y{t)  represents  u{t),  v(t)  or  p(t).  The  wavelet-filtered  time  record  can 
be  reconstructed  by  computing  the  inverse  fast  wavelet  transform  (IFWT) 
of  the  remaining  wavelet  coefficients.  The  discussion  will  focus  on  the  u(t) 
time  record.  However,  all  three  time  records  were  analyzed. 

The  performance  of  the  described  wavelet  filter  depends  on  the  wavelet 
basis  and  the  value  of  k/.  There  are  many  families  of  wavelets  that  could  be 
chosen  for  the  implementation  of  the  wavelet  filter. 

Let  Uu,(t)  be  the  wavelet  filtered  sequence  of  u(t),  then, 

u(t)  =  u^(t)  +  Uu,„(t)  (3.2) 

where  Uu,„(f)  represents  the  noise  present  in  the  filtered  sequence.  As  a 
first  approximation,  assume  that  the  overall  shape  of  the  PSD  of  u(t)  is 
determined  mainly  by  the  presence  of  organized  turbulent  structures.  A 
wavelet  selection  is  made  by  choosing  the  wavelet  basis  for  which  the  PSD  of 
Utu(t)  most  closely  matches  the  PSD  of  u(t).  Three  families  of  wavelets  will 
be  examined. 

The  first  example  uses  a  wavelet  filter  constructed  with  the  8-coefficient 
Daubechies  Extremal  Phase  (D8)  wavelet  basis  (cf.  [8,  page  194]).  The  PSD’s 
of  u(t),  u^{t)  and  Ua,„(t)  for  kf  =  0.9  are  shown  in  figure  3.1. 

For  the  second  example,  the  wavelet  filter  is  implemented  with  the  18- 
coefficient  Coiflet  (C18)  wavelet  basis  (cf.  [8,  page  261])  and  kf  =  0.9.  The 
corresponding  PSD’s  are  shown  in  figure  3.2. 

For  the  third  example,  the  wavelet  filter  is  implemented  with  the  8- 
coefficient  Daubechies  Least  Asymmetric  (L8)  wavelet  basis  (cf.  [8,  page  198]) 
and  kf  =  0.9.  The  corresponding  PSD’s  are  shown  in  figure  3.3. 

The  three  examples  show  that  the  overall  shapes  of  the  PSD’s  of  u{t) 
and  Uy,{t)  are  not  preserved  in  the  high  frequency  range.  However,  it  should 
be  noted  that  the  PSD’s  of  u^{t)  in  the  three  examples  have  similar  shapes. 
It  appears  that  the  PSD  of  the  wavelet-filtered  record  is  independent  of  the 
choice  of  wavelet  basis. 


Figure  3.1:  PSD  of  u(i),  Uy,{t)  and  The  sequence  Uy,{t)  was  obtained 

using  the  D8  wavelet  filter  with  kf  =  0.9 


Figure  3.2:  PSD  of  u{t),  liv/it)  and  The  sequence  u^{t)  was  obtained 

using  the  C18  wavelet  filter  with  kj  =  0.9 
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Figure  3.3:  PSD  of  u{t),  Uy}{t)  and  The  sequence  u^u{t)  was  obtained 

using  the  L8  wavelet  filter  with  kf  =  0.9 


It  is  instructive  to  examine  the  data  compression  achieved  by  the  wavelet 
filter.  From  the  262144  wavelet  coefficients  obtained  from  the  zero  padded 
u{t)  time  record,  only  20373  coefficients  are  used  to  reconstruct  u^{t)  in  the 
first  example,  19672  coefficients  are  used  in  the  second  example,  and  20367  in 
the  third  example.  In  all  cases,  les':  than  10%  of  the  coefficients  are  retained, 
therefore,  significant  data  compression  has  been  achieved.  The  fact  that  the 
number  of  wavelet  coefficients  is  lower  in  the  second  example  is  an  indication 
that  the  signal  energy  for  this  case  is  slightly  more  concentrated  in  wavelet 
space. 

Segments  of  the  u{t)  and  u«,(t)  data  records  are  shown  in  figure  3.4  for  the 
three  examples.  There  is  little  difference  between  the  original  and  the  filtered 
records  because  only  10  percent  of  the  original  energy  has  been  removed. 

These  conventional  wavelet  filters  have  effects  similar  to  that  of  a  low-pass 
filter,  \  alues  of  k/  closer  to  one  cause  a  shift  of  the  equivalent  cutoff  frequency 
towards  higher  frequencies.  Lower  values  of  kf  cause  a  shift  of  the  equivalent 
cutoff  frequency  towards  lower  frequencies.  The  natural  question  that  arises 
is:  what  is  the  difference  between  this  filter  and  the  traditional  low-pass  filters 
in  the  frequency  domain?  One  way  to  resolve  this  issue  is  to  select  a  classical 
filter  that  would  produce  effects  similar  to  the  ones  seen  in  figures  3.1  to  3.3. 
A  lowpass  Butterworth  filter  of  first  order  with  cutoff  frequency  at  2  KHz 
will  produce  the  least  steep  rolloff  in  the  PSD  of  the  filtered  sequence.  The 
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Figure  3.4:  Segments  of  the  u{t)  and  n«,(t)  time  records  for  the  wavelet  filters 
described  in  examples  1  (D8),  2  (C18)  and  3  (L8). 


Figure  3.5:  PSD  of  u(t),  u/(t)  and  u/„(t).  The  sequence  u/(t)  was  obtained 
using  a  first  order  Butterworth  filter  (double  pass  for  phase  compensation) 
with  cutoff  frequency  at  2000  Hz. 
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filter  order  is  effectively  doubled  by  forward  and  reverse  filtering  to  provide 
phase  compensation.  Using  a  notation  consistent  with  (3.2),  let  Uf{t)  be  the 
butterworth  filtered  sequence  of  u(t),  then, 

u{t)  =  Uf{t)  +  Uf^{t)  (3.3) 

where  u/„(t)  represents  the  noise  present  in  the  filtered  sequence.  The  re¬ 
sulting  PSD’s  are  shown  in  figure  3,5. 

Butterworth  filters  of  higher  order  produce  a  rolloff  steeper  than  the  one 
shown  in  figure  3.5.  This  effect  is  not  desired  if  the  shape  of  the  PSD  is  to 
be  preserved.  An  increase  in  the  cutoff  frequency  causes  less  distortion  in 
the  PSD,  however,  the  effect  of  filtering  decreases  because  the  noise  remains 
present  in  the  filtered  record. 

The  wavelet  filter  appears  to  affect  the  data  in  a  range  of  frequencies 
that  vary  from  low  frequencies  to  high  frequencies,  while  the  effect  of  the 
Butterworth  filter  is  concentrated  in  the  higher  frequencies. 


3.2  Weighted  Wavelet  Filtering 

To  extend  the  effects  of  the  wavelet  filter  to  a  wider  frequency  range,  a 
weight  w(m)  can  be  assigned  to  the  rows  of  the  wavelet  coefficients  so  that 
those  coefficients  associated  with  a  common  scaling  factor  a  =  2^  have  the 
same  weight.  For  instance,  in  a  sequence  having  2^  samples,  the  weight 
could  be  set  to  u;(m)  =  2““^.  Figure  3.6  shows  three  examples  of  the  effect 
of  multiplying  the  D8  FWT  coefficients  of  u{t)  by  its  associated  weights 
w{m)  =  1,  w(m)  =  2""*/^  and  2“"*  (using  a  =  0,  a  =  1/2  and  o  =  1, 
respectively).  Similar  results  were  found  using  C18  and  L8  FWT’s. 

The  modified  wavelet  filter  is  constructed  as  it  is  initially  described,  with 
the  exception  that  the  sorted  descending  ordering  is  applied  to  the  weighted 
wavelet  coefficients.  It  should  be  noticed  that  the  weighting  scheme  is  used 
only  to  determine  which  coefficients  are  selected.  The  selected  coefficients  do 
not  undergo  any  modification. 

In  order  to  examine  the  effects  of  the  weighting  scheme  on  the  selection 
of  the  coefficients  Ct,  the  segment  of  the  total  wavelet  power  contained  in 
each  row  m  of  the  coefficient  expansion  will  be  evaluated.  The  weighting 
10(771)  alters  the  selection  of  the  individual  terms  in  each  row  and  thus  the 
frequency  content  of  the  filtered  signal.  Figures  3.7,  3.8  and  3.9  show  the 


24 


0  (LS  1  \S  _  2  3  3^ 

WsMltl  ooiMoiifil  indiK  mIO* 

Figure  3.6:  Effects  of  the  weights  w{m)  =  1  (top),  w{m)  =  (middle) 
and  w{m)  =  2“"*  (bottom)  applied  to  the  D8  FWT  coefficients  of  u{t).  A 
segment  of  the  u{t)  record  containing  7^=  32768  was  used  in  the  computa¬ 
tions. 

signal  energy  versus  m  for  the  three  examples  of  figure  3.6.  The  filter  was 
constructed  for  kf  =  1,  .9.,  .7,  and  ,  .5, 

As  previously  stated,  the  wavelet  coefficients  that  dominate  the  signal 
are  terms  in  rows  2  <  m  <  8  with  center  frequencies  2181  >  /c  >  34, 
respectively.  The  results  indicate  that: 

•  The  weighting  w{Tn)  =  1  (conventional  filtering)  crudely  approximates 
a  low-pass  filter  —  as  previously  discussed  —  by  mostly  eliminating 
coefficients  in  rows  2  <  m  <  6  associated  with  wavelets  that  have  high 
frequency  components. 

•  The  weighting  w{m)  =  2’'”*/^  approximates  a  bandpass  filter  by  elim¬ 
inating  coefficients  associated  with  wavelets  containing  both  high  and 
low  frequency  components. 

•  The  weighting  w{m)  =  2~"*  approximates  a  high-pass  filter  by  elimi¬ 
nating  coefficients  in  rows  5  <  m  <  8  associated  with  wavelets  that 
have  low-frequency  components. 
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Figure  3.7:  Normalized  D8  Wavelet  Power  versus  m  for  w{m) 
and  kj  =  1,  .9,  .7  and  .5 


Figure  3.8:  Normalized  D8  Wavelet  Power  versus  m  for  ii;(m)  = 
1/2)  and  kj  =  1,  .9,  .7  and  .5 
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Figure  3.9:  Normalized  D8  Wavelet  Power  versus  m  for  u;(m)  =  2“”*  (a  =  1) 
and  kf  =  1,  .9,  .7  and  .5 

This  is  more  clearly  displayed  in  figures  3.10,  3.11  and  3.12  where  the 
PSD’s  of  the  reconstructed  filtered  signals  are  shown.  These  figures  cor¬ 
respond  to  the  three  weighting  schemes  for  kf  =  .5.  Similar  results  were 
obtained  for  other  values  of  kf. 

As  a  first  attempt  at  designing  a  wavelet  filter,  extensive  computations 
were  made  utilizing  the  three  weighting  schemes  discussed  above.  The  tech¬ 
niques  were  applied  simultaneously  to  u(t),  v(t),  and  p(t)  time  records  at 
several  locations  across  the  turbulent  boundary  layer. 

A  wavelet  filter  constructed  using  D8  with  w(m)  =  and  kf  =  0.7 
will  be  used  to  illustrate  the  effect  of  weighted  filtering.  This  filter  was 
chosen  because  it  best  preserves  the  overall  spectrum  shape  while  eliminating 
energy  from  both  low  and  high  ends  of  the  spectrum.  The  time  record  of 
unfiltered  and  filtered  signals  for  u(t),  v(t)  and  p(t)  at  the  near  wall  location 
yfb  =  0.08  are  shown  in  figure  3.13.  The  filtered  sequences  are  smoother 
than  the  original  signal,  primarily  because  of  the  low-pass  effect  of  the  filter. 
It  also  appears  that  the  peak  events  are  preserved  in  the  filtered  record.  The 
next  Chapter  will  perform  further  anals^es  of  this  filter  scheme  based  on 
conditional  sampling  of  both  the  filtered  and  unfiltered  time  records. 
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Figure  3.10:  PSD  of  u{t)  (solid)  and  v^(t)  (dashed)  for  k 

(a  =  0) 


Figure  3.11:  PSD  of  u{t)  (solid)  and  u„,{t)  (dashed)  for 

2-m/2  ^  1/2) 
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Figure  3.12:  PSD  of  u(i)  (solid)  and  u^{t)  (dashed)  for  kf  =  .5  and  w(m)  = 
2-”*  (a  =  1) 


Figure  3.13:  Segments  of  u(t),  ittvit),  v(t),  v„,(t),  p(t)  and  Pw(t).  The  se¬ 
quences  u„(t),  v„,(t)  and  Pu,(t)  were  obtained  using  the  D8  modified  wavelet 
filter  with  kf  =  .7  and  ti;(m)  = 
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Chapter  4 

EVALUATION  OF  EVENT 
DETECTION 


An  initial  attempt  to  evaluate  the  filtered  results  will  be  based  on  event 
detection  and  conditional  sampling.  A  preliminary  study  of  the  effect  of  the 
wavelet  filter  on  the  VITA+  slope  event  detection  method  will  be  presented. 
The  VITA  events  of  the  u{t)  sequence  require  the  evaluation  of  (4.1). 


1  o 

u(to,r„)  =  yJ  Zk  “  (*o  +  T)dr 

A  VITA  event  is  detected  if  u{to,T„)  >  The  VITA  parameters 

used  in  this  Chapter  are  =  1,  and  Ty  corresponding  to  24  consecutive 
samples. 

The  +slope  events  are  detected  when  the  following  condition  is  true 

—  I  u(to  +  r)dr  -  ^  w(<o  +  T)dT  >  0  (4.2) 

The  T,  parameter  was  set  equal  to  Ty. 

VITA+slope  events  are  usually  found  by  performing  an  AND  operation 
of  the  VITA  events  with  the  +slope  events.  However,  in  order  to  avoid  the 
double  count  of  a  VITA  event  that  coexists  with  two  different  +slope  events 
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Figure  4.1:  Segment  of  u{t)  (upper),  VITA+slope  events  detected  on  u{t) 
(middle),  and  VITA+slope  events  detected  on  u^{t). 


(or  vice-versa),  a  VITA+slope  event  is  defined  as  a  +slope  event  that  coexists 
with  one  or  more  VITA  events. 

The  application  of  the  VITA+slope  detection  scheme  results  in  the  de¬ 
tection  of  1568  events  when  it  is  applied  to  the  u{t)  sequence,  and  1576 
events  when  it  is  applied  to  the  u^{t)  sequence.  Figure  4.1  shows  a  section  of 
the  u(t)  record  and  the  VITA+slope  events  detected  using  the  original  and 
filtered  sequences. 

Let  q2{t)  and  94 (t)  be  defined  by 


,  ^  _  J  u{t)v{t)  if  u(t)  <  0  and  v{t)  >  0 
^  ~  (  0  otherwise 

,  .  _  f  u{t)v{t)  if  u(t)  >  0  and  v{t)  <  0 
^  ~  I  0  otherwise 


(4.3) 

(4.4) 


The  average  shape  of  the  VITA+slope  events  detected  on  the  u{t)  record 
can  be  found  by  averaging  all  the  events  centered  within  a  window  of  con¬ 
stant  size.  The  average  shapes  of  the  ti(<),  i;(<),  uv{t),  q2{t)  and  g4(t)  ob¬ 
tained  in  this  way  are  plotted  in  figures  4.2  and  4.3.  The  equivalent  pro¬ 
cessing  performed  in  the  filtered  sequences  (using  the  events  detected  by  the 
VITA+slope  scheme  on  the  u^{t)  record)  are  given  in  figures  4.4  and  4.5. 
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Figure  4.2;  Aver^^c  shape  of  p(t),u(t),  t;(t)  and  uv{t)  for  VITA+slope  events 
detected  on  v.^t^ 
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Figure  4.3:  Average  shape  of  p{t),q2{t),  qi{t)  and  uv(t)  for  the  conditional 
sampling  scheme  performed  on  the  VITA+slope  events  detected  on  u{t) 
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Figure  4.4:  Average  shape  of  Pu,(t),Uu„(t),  v^{t)  and  uvv,{t)  =  Uy,{t)v^{t)  for 
VITA+slope  events  detected  on  u^{t). 
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Figure  4.5:  Average  shape  of  Pw{t),q2^{t),  94* (0  and  uv^{t)  =  u^{t)v^{t) 
for  the  conditional  sampling  scheme  performed  on  the  VITA+slope  events 
detected  on  Uv,{t). 


In  all  cases  the  conditional  sampling  results  are  consistent  even  though 
only  70%  of  the  original  signal  energy  is  present  in  the  filtered  records. 
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Chapter  5 

SUMMARY  AND 
DISCUSSION 


The  fundamental  problem  posed  in  this  investigation  is  to  separate  the  signal 
from  the  noise  in  the  temporal  records  defined  in  (1),  (2)  and  (3).  Extensive 
knowledge  of  the  signal  model  is  a  requirement  for  the  design  of  efficient 
filters.  However,  the  signal  model  for  this  problem  is  not  known. 

Conventional  spectral  filtering  does  not  preserve  the  integrity  of  the  signal 
because  of  the  low  signal- to-noise-ratio  in  the  frequency  band  of  interest. 

Some  time-frequency  localization  techniques  have  been  presented  that 
offer  a  distinct  perspective  on  the  formulation  of  the  problem.  The  STFT 
and  the  Morlet,  First  and  Second  Derivatives  of  the  Gaussian  Wavelets  have 
been  applied  to  turbulent  and  wall-pressure  time  records.  The  corresponding 
phase  and  magnitude  spectrograms  and  scalograms  do  not  contain  clearly 
indentifiable  patterns  that  could  lead  to  a  reliable  identification  of  the  signal 
embedd  ?d  in  the  noise. 

The  FWT  is  a  discrete  orthogonal  wavelet  expansion  that  does  not  intro¬ 
duce  the  redundancies  of  both  the  STFT  and  the  continuous  non-orthogonal 
Wavelets.  It  also  has  a  structure  that  can  be  used  for  the  systematic  dis¬ 
crimination  of  frequency  bands. 

A  filtering  tec'.inique  based  on  the  FWT  was  developed.  A  weighting 
scheme  that  allows  a  systematic  selection  of  the  wavelet  coefficients  is  intro¬ 
duced.  One  of  the  effects  of  the  selected  weights  is  to  preserve  the  shape 
of  the  signal  PSD,  which  is  attributed  to  the  existence  of  large  amplitude 
events. 
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An  evaluation  of  the  proposed  filter  was  performed  using  conditional  sam¬ 
pling  based  on  VITA-fslope  detection  on  the  u  filtered  and  unfiltered  records. 
The  results  indicate  that  the  filtered  sequences  preserve  the  events  while 
eliminating  some  of  the  high  and  low  frequency  components.  As  a  result, 
the  process  of  scanning  the  filtered  data  records  becomes  a  possibility  in  the 
search  of  tentative  signal  models. 

Research  is  now  underway  to  examine  other  weighting  techniques  and 
the  subsequent  effect  on  the  filtered  temporal  records.  Furthermore,  in  Part 
II  of  this  research  investigation,  filtering  techniques  will  be  utilized  for  the 
detection  of  groups  or  clusters  of  burst  events  which  arc.  associated  with 
the  passage  of  a  single  large-scale  flow  structure.  This  problem  requires  a 
different  approach  to  wavelet  filtering  and  sampling.  New  data  records  are 
also  being  acquired  for  disturbed  flows. 

ACKNOWLEDGEMENTS 

This  work  was  supported  by  the  Office  of  Naval  Research  with  Dr.  L.  Patrick 
Purtell  as  scientific  officer,  under  Grant  #  N00014-88-K-0141, 

The  data  records  used  in  this  phase  of  the  investigation  were  acquired  by 
Dr.  Vincent  Wilczynski.  The  authors  are  grateful  for  his  cooperation. 


Bibliography 


[1]  Robinson,  S.K.,“  Coherent  Motions  in  the  Turbulent  Boundary  Layer,” 
Ann.  Rev.  of  Fluid  Mech.,  vol.  23,  pp.  601-639,  1991. 

[2]  Farabee,  T.  M.  and  Casarella,  M.  J.,  “Spectral  Features  of  Wall  Pressure 
Fluctuations  Beneath  Ihrbulent  Boundary  Layers,”  Phys.  Fluids  A,  Vol 
3,  No.  10,  pp.  2410-2420, 1991. 

[3]  Kim,  J.,  “On  the  Structure  of  Pressure  Fluctuations  in  Simulated  Tur¬ 
bulent  Channel  Flow,”  J.  Fluid  Mech.,  No.  205,  pp.  421-451,  1989. 

[4]  Wilczynski,  V.,  Casarella,  M.J,  and  Kammeyer,  M.,  “A  Comparison 
of  Data  on  Intermittent  Turbulent  and  Wall  Pressure  Events,”  ASME 
Winter  Meeting,  NCA-Vol  15,  FED- Vol  168,  1993. 

[5]  Snarski,  S.  R.,  and  Lueptow,  R.M.,  “Wall  Pressure  and  Turbulent  Struc¬ 
tures  in  a  TWbulent  Boundary  Layer  on  a  Cylinder  in  Axial  Flow,” 
Submitted  to  J.  Fluid  Mech.,  1994. 

[6]  Morrison,  J.  F.,  Tsai,  H.M.,  Bradshaw,  P.,  “Conditional-sampling 
Schemes  for  Turbulent  flow,  based  on  Variable-interval  time  averaging 
(VITA)  algorithm,”  Experiments  in  Fluids,  Vol.  7,  pp.  173-189, 1989. 

[7]  Wilczynski,  V.,  Organized  Turbxdent  Structures  and  Their  Induced  Wall 
Pressure  Fluctuations,  PhD  Dissertation,  The  Catholic  University  of 
America,  1992. 

[8]  Daubechies,  I.,  Ten  Lectures  on  Wavelets,  Society  for  Industrial  and 
Applied  Mathematics,  1992. 

[9]  Meneveau,  C.,  “Analysis  of  Turbulence  in  the  Orthonormal  Wavelet 
Representation,”  J.  Fluid  Mech.  No.  232,  pp.  469-520, 1991. 


[10]  Farge,  M.  “Wavelet  Transform  and  their  Applications  to  Turbulence,” 
Annu.  Rev  Fluid  Mech.,  pp.  395-457,  1992. 

[11]  Zubair,  L.,  “Studies  in  turbulence  using  Wavelets  Transformations  for 
Data  Compression  and  Scale  Separation,”  Ph.D.  Dissertation,  Yale, 
1993. 

[12]  Lewalle,  J.,  “Wavelet  Analysis  of  Experimental  Data:  Some  Methods 
and  Underlying  Physics,”  25th.  AIAA  Fluid  Dyn.  Conf.,  AIAA  94-2281, 
June  1994. 

[13]  Liandrat,  J.  and  Moret-Bailly,  F.,  “The  Wavelet  Transform:  Some  Ap¬ 
plications  to  Fluid  Dynamics  and  Turbulence,”  Eur.  J.  Mech.  B/Fluids, 
vol.  9,  No.  1,  pp.  111-119, 1990. 

[14]  Press  H.,  et.  al.,  Numerical  Recipes  in  C,  Cambridge  Press,  1992. 

[15]  Wark,  C.  E.  and  Nagib,  H.  M.,  “Experimental  Investigation  of  Coherent 
Structures  in  Turbulent  Boundary  Layers,”  J.  Fluid  Mech.,  No.  230,  pp. 
183-200, 1991. 

[16]  Kaftori,  D.,  Structures  in  the  Turbulent  Boundary  Layer  and  their  Inter¬ 
action  with  Particles,  PhD  Dissertation,  University  of  California,  Santa 
Barbara  ,1993. 

[17]  Kaftori  D.,  et.  al.,  “Funnel-shaped  vortical  structures  in  wall  turbu¬ 
lence,”  Phys.  Fluids,  Vol.  6,  No.  9,  pp.  3035-3050,  1994. 


38 


